print(z <- nsz(14.084,-3.4107,0.0015,1.8832,timepoints))